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ABSTRACT 

A  simplified  form  of  the  three-degrees-of -freedom  nonlinear 
equations  of  the  3000  ton  Surface  Effect  Ship  (SES)  was  lin- 
earized around  an  equilibrium  point  when  a  thruster  failure 
or  thruster  reversal  occured.   Optimal  control  feedback  gains 
were  calculated  with  the  linearized  system  and  applied  to  the 
nonlinear  system.   Evaluation  of  this  design  was  based  upon 
both  comparison  of  the  real-time  human  operator  response  and 
additionally  a  state  feedback  design  based  upon  a  trial  and 
error  procedure. 
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I.   INTRODUCTION 

Over  the  past  years  a  lot  of  attention  has  been  given  to 
the  Surface  Effect  Ship  (SES)  due  to  its  advantages  over  con- 
ventional ship  designs  [Refs.  1,  2,  and  3] 

Recent  investigations  regarding  water  jet  propulsion  SES 
have  been  conducted  on  the  computer  assisted  mode  of  operation 
either  in  a  normal  mode  of  operation  or  during  a  variety  of 
failures  [Ref.  3],   Failures  considered  include,  among  others, 
complete  loss  of  thrust  and  reversal  of  thrust  of  a  nozzle. 

The  purpose  of  this  thesis  was  to  investigate  the  possi- 
bility of  optimal  maneuvering  of  a  3000  ton  Surface  Effect 
Ship  during  the  above  mentioned  emergencies . 

As  discussed  in  [Ref.  3],  turn  maneuvers  for  the  3KSES 
are  the  result  of  turning  moments  generated  by  modulating  the 
waterjet  thrust  vectors.   This  modulation  is  achieved  by  a 
combination  of  thrust  amplitude  modulation  (controlled  by 
varying  the  flow  rates  of  the  four  water jets  to  create  a 
net  differential  thrust,  i.e.,  by  varying  the  power  turbine 
speeds  of  the  individual  propulsion  engines) ,  and  thrust 
vectoring  (controlled  by  varying  the  directions  of  the  four 
waterjet  deflection  sleeves). 

In  the  computer  assist  mode,  the  relative  position  of  the 
helm  determines  the  angular  rate  at  which  the  turn  is  nego- 
tiated.  Full  helm  position  corresponds  to  the  maximum  safe 
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turning  rate  for  a  given  ship  speed.   Calibration  of  the  helm 
input  is  done  within  the  control  computer,  which  then  gen- 
erates commands  to  the  individual  control  effectors  to  produce 
the  commanded  turn  rate . 

In  the  manual  mode,  the  relative  helm  position  generates 
commands  for  the  direction  nozzles  directly.   A  common  position 
command  drives  all  four  deflection  sleeves.   Because  of  that, 
the  rate  of  turn  which  corresponds  to  a  particular  helm  posi- 
tion is  dependent  on  the  hydrodynamic  yawing  moment  which  is 
a  function  of  ship  speed.   The  operator  can  also  generate  a 
thrust  differential  turning  moment  in  the  manual  mode  by 
proper  selection  of  the  propulsion  engine  control  levers. 

To  relieve  the  operator  of  the  work  load  related  to  moni- 
toring and  maintaining  ship  heading  over  extended  periods  of 
time,  a  heading  control  autopilot  is  included  in  the  design 
of  the  heading  control  system.   From  an  operational  point  of 
view,  the  heading  autopilot  is  limited  to  a  "heading  hold" 
function,  and  is  used  only  in  conjunction  with  the  computer 
assist  mode  of  operation.   The  net  result  of  this  operation 
is  that  the  heading  autopilot  replaces  the  operator's  hand  on 
the  helm,  while  the  helm  is  still  the  primary  control  input 
device  for  the  control  of  ship  heading. 

This  is  accomplished  by  direct  feedback  of  the  yaw  rate 
to  a  processor  which  also  gets  an  error  channel  from  the  auto- 
pilot through  a  control  filter  using  proportional  plus  inte- 
gral control  of  the  heading  angle. 
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The  fundamental  optimality  requirement  used  in  [Ref .  23 
is  concerned  with  minimizing  fuel  and  effector  duty  cycle 
while  executing  prescribed  maneuvers.   This  was  done  by 
minimization  of  a  function  of  the  state  variables. 

The  approach  used  in  this  work  was  to  achieve  results  by 

tf    T       T 
minimizing  a  functional  of  the  form  /    (x  Qx  +  U  RlOdt 

"CO    ~   r<«  ~      "•*   ~  «• 

to  compute  the  gains  for  the  state  feedback  control  (see 
Appendix  A) . 

The  resulting  linear  regulator  design  was  tested  on  the 
nonlinear  system  to  compare  the  performance  of  the  craft 
during  a  failure  mode  with  the  performance  of  a  human  operator 
using  a  real  time  simulation  [Ref.  M-].   The  performance  of  the 
design  during  a  failure  was  also  compared  with  a  corrective 
algorithm  for  control  found  by  simple  trial  and  error  pro- 
cedures . 
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II.   EQUATIONS  OF  MOTION  OF  THE  THREE -DEGREE - 
OF  FREEDOM  (3DOF)  MODEL 


The  three  degrees  of  freedom  to  be  used  in  this  analysis 
are:   surge,  sway  and  yaw  [Ref.  5].   Therefore,  the  following 
assumptions  and  approximations  were  made: 

1.  Roll  is  zero. 

2.  Pitch  and  draft  remain  constant. 

3.  The  origin  of  the  ship's  coordinate  system  is  located 
at  the  center  of  gravity  of  the  ship. 

4- .   At  time  =  0,  the  ship's  coordinate  system  and  the 

fixed  (geographical)  coordinate  system  have  the  same  origin 

and  the  same  coordinate  axes,  i.e.,  xrTAW  =  x   =0  and  y.TA„  = 

'      '   NAV     o         -'NAV 

y   =  0  (see  Fig.  1)  . 
J  o 

5.  At  time  =  0  the  ship  is  heading  in  the  positive  xMAV 
direction. 

6 .  Drift  angle  subtraction  effect  on  thrust  angle  is  zero 
Then,  as  a  good  approximation,  the  equations  of  motion 

are : 

Surge :   m(u-vr)  =  X  (1) 

Sway:    m(v+ur)  =  Y  (2) 

Yaw:     I  r=N  (3) 
z 

Navigation: 

x   =  ucos^-vsini/;  (M-) 

y   =  usintjj+vcos^  (5) 


13 


'NAV 


x 
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SHIP'S  PATH 
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NAV 


Fig,  1  -  Definition  of  Coordinate  System 

NOTE:   Rudder  condition  for  right  turn  is  shown.   All  runs 
start  at  geographic  origin. 
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V   =  /  u2+v2  (6) 

s 

6  =  tan_1(^)  (7) 

In  the  above  equations: 

m  =  mass  of  the  rigid  ship 

I   =  moment  of  inertia  about  the  z  axis 
z 

u  =  V  cos3  =  velocity  in  the  x  direction  (Surge) 

v  =  V  sing  =  velocity  in  the  y  direction  (Sway) 

r  =  ty    =  angular  velocity  about  the  z  axis 

ip  =  heading  (yaw)  angle,  Fig.  1 

3  =  drift  angle, Fig .  1 

6  =  rudder  angle,  +6  produces  a  left  turn 

Also: 

x      =    /x   dt  y      =    /v   dt 

o  o  Jo  Jo 

u      =    /udt  v      =    /vdt 

r      =    /rdt  ip      =    /ipdt 

The  right-hand  side  of  the  equations  of  motion  have  been 
represented  by  X,  Y,  N.   These  symbols  represent:   summation 
of  forces  (X,  Y)  and  moments  (N) .   Each  summation  consists  of 
many  terms.   For  the  initial  simplified  model,  only  the  few 
terms  that  are  major  contributions  have  been  retained.   Those 
terms  are  considered  sufficient  to  produce  a  simulation  which 
would  be  in  close  agreement  with  obtained  data.   At  a  later 
time,  the  model  may  be  improved  by  adding  more  terms  and  ad- 
justing coefficients. 
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After  reduction  of  the  right-hand-sides,  these  force  and 
moment  summations  become: 

X  =  GTcos6-CDXu2  (8) 

Y  =  GTsinS-CDYv|v|  (9) 

N  =  -(GTsin6)l  +C_vv|vll  (10) 

o   DY  '  '  w 

where  (see  Fig.  2): 

GT  =  gross  thrust  in  lbs. 

6 =6  =5   =  angle  of  all  thrust  vectors  (starboard  and  port) 
s   p      to  r 

in  radians 

Cny  =  coefficient  of  drag  in  the  direction  of  surge 

Cny  =  coefficient  of  drag  in  the  direction  of  sway 

1   =  thrust  level  arm  in  feet 
o 

1   =  sway-drag  level  arm  in  feet 
w      J  & 

Also,  note  that  the  drag  forces  have  been  chosen  as  lumped 
parameters  acting  at  the  two  points  P  £  P' ,  i.e.   P£P'  are 
equivalent  force  centroids . 

The  length  1   is  the  distance  from  the  center  of  gravity  to 

the  stern,  and  1   is  the  distance  (in  the  x-direction)  from  the 

'       w 

center  of  gravity  to  the  points  P  £  P',   1   is  taken  as  a  fixed 

parameter  (it  could  be  charged  by  various  conditions  of  ship 

loading)  but  1   is  adjustable  in  Order  to  be  used  as  a  para- 
w 

meter  to  make  the  model  fit  test  data. 

The  coefficient  C~y  in  eq.  8  is  a  lumped  parameter  equiva- 


lent for  the  hull  friction  effect  (slender  body  theory).   The 


2 
term  C_yu  models  the  hull  drag,  aerodynamic  and  hydrodynamic 
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effects.   The  bubble  drag  is  not  included  in  the  initial 
model. 

In  like  manner,  the  coefficient  C    in  eqs .  (9)  and  (10) 
models  the  cross  flow  drag  effects.   Note  also  that  the  added 
mass  term  has  been  temporarily  omitted  and  will  be  added 
later,  especially  in  the  3KSES  study  since  its  contribution  to 
the  XR-3  model  is  almost  nonexistent  [Ref.  6], 

Substitution  of  equations  (8),  (9)  and  (10)  into  the 

equations  (1),  (2)  and  (3)  provides  the  basic  equations  of 

motion  for  the  three-degrees-of -freedom  model: 

u  =  -  (GTcosS  -  C^vu2)  +  vr  (11) 

m  DX 

v  =  -  (GTsin6  -  C^vlv!)  -  ur  (12) 

m  DY  '  ' 

r  =  y~    (-(GTsin<5)lo  +  CDyv|v|lw)  (13) 

z 

Equations  (11),  (12)  and  (13)  are  applicable  to  any  SES 
craft.   Only  the  particular  parameter  values  to  be  used  in 
the  equations  would  depend  on  the  particular  craft  (XR-3  or 
3KSES)  . 

In  the  case  of  the  3KSES ,  the  equations  take  the  form 

(see  Fig .  3 ) : 

u    =    -    (-CT,vu2+(T7+T8+T9+T10)cos6  )    +    vr 
m  DX 

v    =    -    ((T7+T8+T9+T10)sin6-CTWv|v|  ) -ur 

m  DY    '     ' 

r   =  i-   (CDY(ww)v|v|    +    [T7cos5S4+T8cos5S3 
z 

-   T9cos<5S      -   T10cos5S2    -    (T7sin<5+T8sin6    +   T9sin5 

+   T10sin6 )(00)]    +  A00uv(ww)) 
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\\>    -    r 

x      =   ucosiii-vsin^ 

O  r  r 

y   =  usinijj  +  vcosij; 

For  the  purpose  of  this  study  the  thruster fs  angle  6  was 
the  same  for  all  thrusters  except  the  starboard  one,  which 
in  the  case  of  thrust  failure  was  considered  0°  (thruster 
oriented  straight  ahead)  and  in  the  case  of  thrust  reversal 
it  was  considered  as  being  18  0°  (thruster  oriented  straight 
astern) . 

In  the  first  case,  the  thrust  of  the  failed  thruster  was 
set  to  zero  pounds  and  in  the  second  case  the  thrust  was  set 
to  full  value  (98000  pounds  for  56  Kts  operation  speed).   That 
means  that  either  in  the  case  of  thrust  failure  or  thrust 
reversal  of  the  starboard  thruster,  the  rest  of  the  thruaters 
are  ganged  together  and  act  in  the  same  direction. 

Therefore,  the  angle  5  can  be  redefined  as  Z7,  Z8,  Z9  and 
Z10  as  shown  in  Fig.  3  to  correspond  to  the  appropriate  thrus- 
ters and  in  both  cases  Z7-Z8-=Z9  with  Z10  being  either  0  or  tt 
radians  . 

In  order  to  linearize  the  equations,  the  following  assump- 
tions were  made: 

1.   The  maximum  rudder  angle  is  less  than  or  equal  to  3  0 
degrees  and  for  most  high  speed  manueuvers  the  rudder  angle 
stays  below  15  degrees.   Thus,  by  letting  cos<5~l  and  sin6~6 
(where  6  is  in  radians),  the  difference  from  the  actual  values 
for  <5  =  ,2618  rads  =  15°  is 
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a.  cosl5°    =    .966,    difference    =    1":^6   x   100%=3,5% 

.96  6 

b,  sinl5°  =  .2588,  difference  =  — 6 


.2588 

x  100%  =  1.15% 
Differences  with  these  magnitudes  should  not  introduce 
large  deviations  from  the  linear  design  when  applied  to  the 

nonlinear  equations.   Thus  in  the  analysis  5     was  restricted 

^  J  max 

to  15  degrees  or  if  6     exceeds  15  degrees  it  probably  will 
to  max  &        tr  j 

be  of  short  duration  and  will  then  not  introduce  large  errors 
in  the  design  values  for  control. 

2.   The  yaw  angle  \\>   was  considered  to  be  less  than  10 
degrees  which  agrees  very  well  with  the  results  obtained  from 
the  simulation  of  the  craft. 

3  .   Even  though  the  sign  of  the  sway  velocity  changes  de- 
pending on  the  particular  situation  occuring ,  the  term  v|v[ 

•a  A  2 

was  considered  as  v  . 

With  the  above  assumptions,  the  equations  for  the  3KSES 

become: 

f,  =  u  =  -  (-CTWu2+T7+T8+T9+T10)  +vr 

1  m    DX 

f0  =  v  Z   -  C(T7+T8+T9+T10)6-C^„v2)  -  ur 

2  m  DY 

f3  =  r  =  j-   (CDYv2(ww)+T7S1++T8S3-T9S1 
z 

-  T10S    -(T7+T8+T9+T10)6 (00 ) +A22uv(ww) ) 


fr    =   >:      ~   u-vib 

o  o 
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f  _  =  y   ~  ml;  +  v 
6    ^o     r 


Defining  the  state  variables 


x, -u-u 
1     o 


x2=v-vQ 


x»=r-r 
3    o 


X4=^o 


Xr=X  -X 

5   o   oo 


x  =y  -y 
6  Jo    Joo 


M  =T7-T7 
1       o 


U2=Z7-Z7o 


Then,     x,  =  u 


x^  =  v 


x^  =  r 


x, 
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xr  =  x 


xc  =  y 
6    v  o 


The  linear  state  equation  coefficients  are 
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In  standard  form  the  system  becomes 

x  =  Ax  +  B  P  where 

~    «>•  «>•    ~  x  ~ 
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So,    by    linearizing   around   the    equilibrium   point    (u    ,    v    , 

r    ,    4>    ,    x       ,    y      )    the    equations    become: 
o        o        oo      7oo  ^ 

-2Cnvu  . 

DA    O  ,      1    n 

xn     =    x,  +    r   x„    +    v   x_    +    —  C. 

1  m  1  o2  o3ml 

2Cnvv  Z7 
a                                     DY   o  o    n 

x~    =    -r   x,     -   x0    -   u  x,    +  C-, 

2  ol  m  2         •  o    3  ml 


T7    +T8+T9 

+  (— ^ — >  c? 

m  2 

A00v    (ww)                 (2Cnvv    +A00u    )(ww) 
11    Q ,  DY    O       2  2    O 

X3  I  Xl  I  2 


S, -Z7    (00)  (00)(T7    +T8+T9) 

♦  (JUjfi ,  Cl j2 c2 


X4    =    X3 


xc     =     X.      -    ll)    x_     -    V    X. 

5  1         To    2  o    4 


X6    =    Vl    +    x2    +    uoX4 


28 


The  above  equations  were  used  to  determine  by  computer 
simulation  the  required  gains  for  the  state  feedbacks  using 
the  linear  regulator  solution  [Refs.  7,  8  and  9]. 

In  order  to  obtain  the  equilibrium  point,  the  nonlinear 
equations  were  used  by  setting  u=0,  v  =  0,  r=0  and  determining, 
after  a  failure  had  occured  and  the  system  being  corrected  as 
far  as  was  possible  by  the  trial  and  error  algorithm  (see 
Chapter  III),  the  "steady  state"  values  for  the  above  men- 
tioned variables. 

Two  "steady  state"  runs  were  made,  one  with  the  starboard 
thruster  failed  and  the  other  with  the  starboard  thruster 
reversed . 

The  "steady  state"  values  obtained  are  listed  in  Tables 
I  and  II. 

TABLE  I 
Steady  State  Values  (Thruster  Failure) 

Surge  Velocity  Sway  Velocity  Angular  Velocity  Yaw  Angle  Track 

Deviation 

(rad/s)  (rad)      (ft) 

-.7415xl0~7       -.0254     10.88 


(ft/s) 

(ft/s) 

7^.422 

1.86 

Port  Thrust 

(lbf) 

87600 
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TABLE  II 

Steady  State  Values  (Thruster  Reversal) 

Surge  Velocity  Sway  Velocity  Angular  Velocity  Yaw  Angle  Track 

Deviation 


(ft/s) 

(ft/s) 

(rad/s) 

(rad) 

(ft) 

55.66 

2.62 

-.3417xl0~5 

-.04709 

21.41 

For  this  study,  the  steady  state  value  achieved  by  the 
simulation  with  the  thrust  failure  were  used  for  the  develop- 
ment of  the  linear  system  and  the  subsequent  use  of  the  opti- 
mal control  scheme.   The  choice  was  quite  arbitrary  but, 
nevertheless,  the  results  show  that  the  system  does  behave 
well,  no  matter  which  one  of  the  two  emergencies  (thrust 
failure  or  thrust  reversal)  occurs. 
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III.   TRIAL  AND  ERROR  DESIGN 

A.  INTRODUCTION 

Taking  into  account  the  impressive  speed  of  the  surface 
effect  ship,  it  is  evident  that  there  is  a  need  for  controlling 
it  either  during  straight  runs  or  during  an  emergency  which 
would  be  manifested  by  the  loss  of  thrust  of  one  thruster  or 
even  by  the  thrust  reversal. 

As  was  seen  from  the  results  of  computer  runs  with  the 
model,  the  behavior  of  the  ship  during  a  failure  was  pre- 
dictable but  nevertheless  disturbing  because  the  ship  veered 
quite  rapidly  from  the  track  that  it  was  following  before 
the  failure. 

Therefore,  the  goal  of  the  control  procedure  was  twofold. 
First,  it  was  required  to  maintain  control  of  the  ship  and 
second  to  maintain  that  control  in  some  acceptable  manner. 

B.  ANALYSIS 

As  was  evidenced  (see  Fig.  4),  when   the  outer  starboard 
thruster  experienced  a  complete  thruster  failure,  the  ship, 
without  any  corrective  action  with  respect  to  nozzle  movement 
or  thrust  variation  of  the  remaining  jet  nozzles,  veered  to 
the  right  and  actually  entered  a  turn. 

By  analyzing  this  behavior,  it  was  seen  that  the  para- 
meters changing  the  most  were  the  yaw  angle  ip  and  the  deflec- 
tion from  the  original  track  which  was  assumed  straight  ahead. 
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Intuitively  then,  it  was  evident  that  from  the  available 
variables,  ty   and  y   would  be  most  eligible  to  control  in  order 
to  bring  the  ship  back  or  close  to  the  original  track.   Be- 
cause it  was  quite  clear  that,  if  the  heading  (yaw  angle)  was 
minimized  as  well  as  the  deviation  from  the  track,  the  ship 
would  follow  again  the  original  track  or  a  track  that  was 
parallel  to  that  one. 

Therefore,  the  corrective  action  chosen  was  nozzle  de- 
flection as  a  function  of  the  yaw  angle  and  the  deviation 
from  the  original  track. 

Furthermore,  since  positive  nozzle  angle  drives  the  ship 
to  the  left,  by  selecting  deviation  to  the  right  as  being 
positive  and  yaw  angle  to  the  right  being  positive  to  the 
right,  it  was  quite  clear  that  the  nozzle  action  would  only 
involve  ip  and  y   without  any  sign  change. 

Thus ,  the  problem  reduced  to  the  point  of  determining  the 
appropriate  gains  for  ip  and  y  and  the  rudder  action  took  the 
form  z=K-^+K9y  . 

C.   TRIAL  AND  ERROR  PROCEDURE 

As  a  beginning,  K,  was  arbitrarily  selected  to  be  1000 
and  K_  being  5 ,   It  was  seen  that  the  action  of  the  rudder 
was  too  heavy  and  the  movement  of  the  ship  was  quite  erratic 
and  at  the  end  became  unbounded . 

Therefore,  both  gains  were  decreased  to  K  =8  and  K  =,035, 
which  gave  a  more  acceptable  performance  (Fig,  14). 
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The  thus  established  rudder  algorithm  proved  to  be  adequate 
even  when  a  reversal  failure  of  the  starboard  thruster  was 
inflicted  (Fig.  21). 

As  a  further  improvement,  the  idea  of  modulating  the 
thrust  of  the  left  thruster  during  a  failure  was  brought  into 
effect.   That  is,  the  thrust  was  varied  to  smaller  values  by 
the  amount  of  the  deflection  of  the  ship  off  the  original 
track  (Fig.  20). 

The  above  scheme  was  used  for  all  the  runs  presented  in 
Figs.  14  through  27. 

In  order  to  evaluate  all  the  above  corrective  actions, 
the  performance  indices  described  in  Chapter  IV  were  used. 

As  can  be  seen,  the  ship  actually  moves  down  the  track 
with  a  large  drift  angle.   That  is,  although  a  corrective 
rudder  action  is  maintained,  the  ship  does  not  turn  into  a 
circle  and  also  it  does  not  veer  to  the  left  or  to  the  right 
although  there  exists  a  sway  velocity  and  asymmetrical  thrust. 
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IV.   LINEAR  REGULATOR  DESIGN 

A.  INTRODUCTION 

The  theory  behind  the  linear  regulator  is  well  established 
and  need  not  be  repeated  here.   For  a  brief  discussion,  see 
Appendix  A. 

After  the  linearization  of  the  system's  equations  around 
an  equilibrium  point,  the  IBM/ 3 60  was  used  to  obtain  the 
solution  of  the  matrix  Riccati  equation,  using  the  method 
presented  in  [Ref s .  8  and  9]. 

B.  ANALYSIS 

It  was  evident  that  in  order  to  achieve  a  parallel  track 
for  the  craft  after  a  failure  had  occured,  the  deflection 
from  the  original  track  should  come  into  play.   Also,  in 
order  for  the  craft  to  maintain  this  parallel  track,  the  yaw 
angle  should  be  minimized  as  well  as  the  angular  velocity 
about  the  z  axis. 

As  with  the  trial  and  error  procedure  the  corrective 
action  chosen  was  ganged  deflection  of  the  operating  nozzles 
and  thrust  control  of  the  port  thruster  to  compensate  for  the 
failed  starboard  thruster. 

C.  PERFORMANCE  INDICES 

One  of  the  basic  problems  of  modern  control  theory  is 
the  translation  of  system  specifications ,  often  in  such 
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subjective  terms  as  "good  rise  time  with  reasonable  overshoot", 
into  the  form  of  a  performance  index.   Furthermore,  after  the 
selection  of  a  performance  index,  the  designer  has  to  make 
an  initial  guess  for  the  weighting  of  the  state  variables  and 
adjust  them  after  comparing  time  histories  of  the  controlled 
system. 

For  the  purpose  of  this  study,  the  following  performance 
index  was  selected: 


J  =  /tf  (xTQx  +  pTRp)  dt 


Concerning  the  weighting  of  the  variables,  considerable  ex- 
perimentation was  done.   In  the  final  form  of  the  above 
performance  index  the  weighting  matrices  Q  and  R  were: 


Q  = 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

20000 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.035 

R  = 
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The  final  weights  used  reflect  the  importance  given  to  the 
yaw  angle,  the  angular  velocity  about  the  z  axis  and  the 
craft  deviation  from  the  original  track. 

The  optimal  control  inputs  obtained  are: 


0 


-.0316 


0 


.681+ 


-.00046       -.0007         0 


185.7 


110.7 


0 


.186 


The  response  of  the  nonlinear  system  with  the  above 
feedback  is  shown  in  Figs.  2  8  through  41, 

In  order  to  compare  the  performance  of  the  craft  under 
the  trial  and  error  algorithm,  the  optimal  algorithm  and  the 
real  time  simulation,  the  following  two  additional  perfor- 
mance indices  were  used: 


J,  =  /tf  y  2dt   and   J0  =  /tf  (Z7)2dt. 
1        o  l 
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V.   COMPARISON  OF  RESULTS 

A.   THRUST  FAILURE  MODE  (OPTIMAL  ALGORITHM  VS.  TRIAL 
AND  ERROR  ALGORITHM) 

Observation  of  the  sea  track  (x   vs .  y  )  (Figs.  14  and  28) 

o      o     ° 

indicates  that  the  trial  and  error  algorithm  is  much  less  of 
the  desired  control  algorithm  when  compared  to  the  optimal, 
because  it  appears  that  the  trial  and  error  algorithm  does 
not  react  fast  enough  to  stabilize  the  vehicle  whereas  the 
optimal  does  so  immediately.   The  trial  and  error  algorithm 
does  seek  out  the  steady  state  value  of  the  track  deflection 
but  the  comparison  of  the  two  tracks  is  dramatic  in  showing 
how  much  better  off  is  the  optimal  algorithm  in  controlling 
the  craft  during  a  thrust  failure.   It  must  be  kept  in 
mind  in  making  this  observation  that  the  trial  and  error 
design  used  in  this  comparison  was  not  developed  to  its  full 
capabilities.   As  stated  earlier,  the  trial  and  error  design 
was  intended  only  as  a  means  of  establishing  the  linear 
system  design. 

From  the  time  history  shown  in  Figs,  1M-  and  28,  the  vehicle 
for  the  same  failure  scenario  reaches  a  steady  state  point  with 
the  optimal  control  in  22  seconds  whereas  with  the  trial  and 
error  control  the  steady  state  point  would  be  reached  some- 
where around  2  00  seconds  or  even  greater,  due  to  the  lightly 
damped  control  system. 
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As  can  be  seen  from  Figs.  2  0  and  3  5  the  port  thruster  is 
much  more  active  with  the  trial  and  error  algorithm.   The 
range  of  modulation  covers  about  14000  pounds  whereas  utiliza- 
tion of  the  optimal  algorithm  gives  a  total  thrust  modulation 
of  2  pounds.   This  does  not  reflect  exactly  what  happens  to 
the  port  thruster,  because  as  stated  earlier  in  Chapter  II  as 
soon  as  the  failure  occurs,  the  port  thrust  drops  to  90%  of 
its  full  value. 

It  should  also  be  noted  that  the  trial  and  error  solution 
seeks  the  steady  state  value  of  the  port  thruster  obtained  by 
the  optimal  solution. 

Comparing  Figs.  18  and  32  it  can  be  seen  in  Fig.  32  that 
the  forward  and  side  thrust  as  well  as  the  yaw  moment  reach 
a  steady  state  within  10  seconds  after  the  initiation  of  the 
thruster  failure  when  the  optimal  solution  is  used.   Instead, 
when  the  trial  and  error  solution  is  used,  the  craft  as  shown 
in  Fig.  18  seeks  the  same  steady  state  but  even  at  100  sec 
after  the  initiation  of  the  failure  it  has  not  stabilized.   On 
the  average ,  the  values  obtained  from  both  algorithms  are 
approximately  the  same. 

However,  it  is  readily  apparent  that  the  trial  and  error 
control  system  has  poles  close  to  the  joy-axis  and  as  a  result 
exhibits  a  response  of  a  lightly  damped  system  whereas  the 
optimal  control  system  responds  with  much  faster  time  constants, 
which  would  indicate  that  this  system  has  poles  at  a  greater 
distance  from  the  jco-axis. 
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B.   THRUST  REVERSAL  MODE  (OPTIMAL  ALGORITHM  VS.  TRIAL 
AND  ERROR  ALGORITHM) 

All  the  previous  comments  that  were  made  for  the  thrust 
failure  mode  also  apply  to  the  case  of  a  thrust  reversal  for 
the  starboard  thruster.   However,  due  to  the  increased 
severity  of  the  forces  acting  during  the  thrust  reversal  and  the 
greater  magnitude  of  the  yaw  moment  (see  Figs,  3  2  and  39), 
the  optimal  control  exhibits  a  slightly  longer  settling  time 
(Figs .  28  and  35) . 

During  the  thrust  reversal  scenario,  the  absolute  magni- 
tudes of  the  state  variables  is  larger  than  those  experienced 
during  the  thrust  failure  case. 

Comparison  between  the  optimal  solution  time  histories 
and  trial  and  error  solution  time  histories  again  confirms 
the  fact  that  optimal  state  variable  feedback  provides  a  more 
highly  damped  system,  with  steady  state  values  being  reached 
within  20  seconds  from  the  initiation  of  the  failure  (see 
Figs.  36,  38,  39  and  40).   The  trial  and  error  solution  is 
shown  in  Figs.  22,  24,  25  and  26  which  again  exhibits  a 
lightly  damped  time  history  which  would  not  be  acceptable 
either  from  a  mechanical  or  ship's  efficiency  viewpoint. 

It  is  observed  (Fig,  23)  that  the  trial  and  error  solution 
requires  the  operational  thrusters  to  be  extended  to  their 
mechanical  limits  of  + . 5  radians  even  at  the  time  =100  seconds 
point,  whereas  the  optimal  solution  in  Fig,  37  has  a  peak 
value  for  the  thrusters  of  approximately  +.4  radians  and 
never  requires  them  to  extend  to  their  mechanical  limit. 
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Fig,  4  -  Sea  Track,  Thruster  Failure,  No  Correction 
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Fig.  5  -  Surge  and  Sway  Velocities  Vs.  Time,  Thruster  Failure, 
No  Correction 
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Fig.  6  -  Yaw  Rate,  Heading  and  Slip  Angle  Vs.  Time,  Thruster 
Failure,  No  Correction 
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Fig.  7  -  Forward  Force,  Side  Force  and  Turn  Moment  Vs.  Time, 
Thruster  Failure,  No  Correction 
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Fig.  8  -  Surge,  Sway  and  Yaw  Acceleration  Vs.  Time,  Thruster 
Failure,  No  Correction 
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Fig,  9  -  Sea  Track,  Thruster  Reversal,  No  Correction 
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Fig,  10  -  Surge  and  Sway  Velocities  Vs.  Time,  Thruster 
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Fig.  11  -  Yaw  Rate,  Heading  and  Slip  Angle  Vs.  Time,  Thruster 
Reversal,  No  Correction 
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Fig.  12  -  Forward  Force,  Side  Force  and  Turn  Moment  Vs.  Time, 
Thruster  Reversal,  No  Correction 
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Fig,  13  -  Surge,  Sway  and  Yaw  Acceleration,  Thruster  Reversal, 
No  Correction 
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g,  14  -  Sea  Track,  Thruster  Failure,  Trial  and  Error  Solution 
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Fig.  15  -  Surge  and  Sway  Velocities  Vs.  Time,  Thruster 
Failure,  Trial  and  Error  Solution 
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Fig.  16  -  Effector  Action  Vs,  Time,  Thruster  Failure,  Trial 
and  Error  Solution 
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Fig.  17  -  Yaw  Rate,  Heading  and  Slip  Angle  Vs.  Time,  Thruster 
Failure ,  Trial  and  Error  Solution 


53 


o 

- 


0C"000e        00-09S2       00'0252        CO'0882       00"0ii82        00"0C82        00*09£2        00"02i2        GO'Gesf 


(jqi)  OIx  soj:oj  papwaoj 


0CT02I  00*501 


00'OS 


00'SZ. 


CO'OS 


OCTSfi 


3C  *0£ 


oo  •  =  : 


00*0 


(jqi)80ix  aojtoj   apxs 


CG"08 


Fig.    II 


)C-G< 


OP  ")h 


00'02  00*0-  00*02- 

"n 


30'0n-         00*09-         00"08- 


-  Forward  Force,  Side  Force  and  Turn  Moment  Vs.  Time, 
Thruster  Failure,  Trial  and  Error  Solution 


54 


80  -0 


os-s 


OCTO  80*0- 


GC'S 


01*0  [CO 


91  -0- 


ti2  *0 


2£"0-  Ofi'O- 


(^03S/^J)^l 


3S"I 


DO- 1 


0S*0 


00*0  OS'O- 


00*  I- 


(     OaS/*].J)A 


eo*c  oo'o        20*0- 

(   oa's/p-ea)      oixJ 


os*:- 


SO'O-  £0*0-  01*0- 


Fig .  19  -  Surge,  Sway  and  Yaw  Acceleration  Vs.  Time,  Thruster 
Failure ,  Trial  and  Error  Solution 
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Fig,  20  -  Port  Thruster  Action  vs.  Time,  Thruster  Failure, 
Trial  and  Error  Solution 
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Fig.  21  -  Sea  Track,  Thruster  Reversal,  Trial  and  Error 
Solution 
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Reversal,  Trial  and  Error  Solution 
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Fig.  23  -  Effector  Action  Vs.  Time,  Thruster  Reversal,  Trial 
and  Error  Solution 
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Fig.  24  -  Yaw  Rate,  Heading  and  Slip  Angle  Vs.  Time,  Thruster 
Reversal ,  Trial  and  Error  Solution 
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Fig.  25  -  Forward  Force,  Side  Force  and  Turn  Moment  Vs.  Time, 
Thruster  Reversal,  Trial  and  Error  Solution 
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Fig.  26  -  Surge,  Sway  and  Yaw  Accelerations  Vs.  Time,  Thruster 
Reversal,  Trial  and  Error  Solution 
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Fig.  27  -  Port  Thruster  Action  Vs.  Time,  Thruster  Reversal, 
Trial  and  Error  Solution 
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Fig.  28  -  Sea  Track,  Thruster  Failure,  Optimal  Solution 
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Fig.  29  -  Surge  and  Sway  Velocities  Vs.  Time,  Thruster 
Failure,  Optimal  Solution 


65 


o 


- 
o 


i-O 


G9*0 


STi 


O 
H 
tSI 


oe-o 


s:  -o 


■y- 

00*0- 


.0 


s 

!  o 

.o 

O 

(1) 

o  w 

=•  s 

H 

Eh 

o 
_d 


r? 


I'O-     C£'0- 


sfi'o-   os-o- 


(P^UZ 


02"  I 


06' 


CS-Q 


0£-Q 


0£"0- 

(ppa)QIZ 


:s-!- 


Fig .  30  -  Effector  Action  Vs.  Time,  Thruster  Failure,  Optimal 
Solution 
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Fig.  3  2  -  Forward  Force,  Side  Force  and  Turn  Moment  Vs.  Time, 
Thruster  Failure,  Optimal  Solution 
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Fig.  34  -  Port  Thruster  Action  Vs.  Time,  Thruster  Failure, 
Optimal  Solution 
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Fig.  35  -  Sea  Track,  Thruster  Reversal,  Optimal  Solution 
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Fig.  36  -  Surge  and  Sway  Velocities  Vs.  Time,  Thruster  Reversal, 
Optimal  Solution 
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Fig.  37  -  Effector  Action  Vs,  Time,  Thruster  Reversal, 
Optimal  Solution 
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Fig .  3  9  -  Forward  Force,  Side  Force  and  Turn  Moment  Vs.  Time, 
Thruster  Reversal,  Optimal  Solution 
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Fig.  40  -  Surge,  Sway  and  Yaw  Accelerations  Vs.  Time,  Thruster 
Reversal,  Optimal  Solution 
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Fig.  41  -  Port  Thruster  Action  Vs.  Time,  Thruster  Reversal, 
Optimal  Solution 
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The  sea  track  (see  Fig.  35)  of  the  optimal  solution 
becomes  stabilized  at  approximately  18  feet  to  the  right  of 
the  original  track  (due  north) ,  with  a  slight  tendency  to 
drift  to  the  right.   This  is  not  of  significant  consideration 
due  to  the  relatively  small  magnitude  of  drift  and  the  fact 
that  during  the  conduct  of  an  actual  emergency  situation 
proper  corrective  action  would  be  also  taken  by  the  crew  of 
the  ship  in  the  first  10  to  15  seconds  of  the  emergency. 

It  is  important  to  keep  in  mind  that  the  ship's  response 
for  both  the  trial  and  error  and  the  optimal  solutions  was 
the  result  of  totally  automatic  operation  and  any  human 
operator  inputs  were  not  considered. 

The  deviation  to  the  right  obtained  by  the  use  of  the  trial 
and  error  algorithm  at  100  seconds  as  shown  in  Fig.  21  has  an 
average  value  of  about  18  feet  and  is  fluctuating  around  this 
value  by  a  magnitude  of  6  feet  which  would  be  more  than  enough 
to  discrupt  normal  ship  activity. 

In  addition,  the  situation  is  more  aggravated  when  surge, 
sway  and  yaw  accelerations  (see  Figs.  2  6  and  33)  are  taken 
into  account  because  even  though  they  are  small,  they  are 
still  significant  after  100  seconds. 
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C.   COMPARISON  OF  THE  OPTIMAL  SOLUTION  AND  THE 
REAL=TIME  HUMAN  OPERATOR  RESPONSE 

As  stated  in  the  introduction ,  a  real  time  simulation 
of  f ive-degree-of-freedom  Surface  Effect  Ship  was  available 
[Ref.  4],  utilizing  a  human  operator  to  input  only  ganged 
thruster  action,  the  thrust  being  constant. 

Comparison  of  the  optimal  solution  and  the  above  men- 
tioned simulation  gave  the  following  results  during  a  thrust 
reversal  scenario: 

1.   Initial  recognition  of  the  emergency  by  the  human 
operator  did  not  occur  until  approximately  three  seconds  into 
the  emergency,  as  .compared  to  the  optimal  solution  with  an 
instantaneous  recognition  of  the  failure.   Consequently, 
as  shown  in  Table  IV,  the  absolute  values  of  all  the  outputs 
reached  by  a  human  operator  controlled  SES  are  much  greater 
than  those  observed  in  the  optimal  control,  given  the  same 
initial  conditions . 

TABLE  III 

Comparison  of  Optimal  Solution  and 
Human  Operator  Response 

HUMAN  OPERATOR  OPTIMAL  SOLUTION 

Time  for  first  overshoot(sec) 

Deflection  from  original 
track  (ft) 

Time  to  settle  back(sec) 

Final  def lection(f t) 

Turn  rate(rad/sec) 


39 

22 

100 

18 

(worst  case) 

78 

35 

0 

18 

.007 

.0065 
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Specific  examples  are  as  follows: 

1.  Maximum  track  deviation  obtained  by  the  ship  with  the 
optimal  control  in  effect  was  approximately  18  feet,  whereas 
with  the  human  operator  maximum  track  deviation  varied  from 
35  feet  to  100  feet,  depending  upon  operator  response  and 
familiarization  with  the  failure  as  well  as  the  ability  of  the 
operator  to  recognize  the  onset  of  an  emergency  situation 
through  the  visual  cues  provided  by  the  real  time  simulation. 

2.  Turn  rate  was  approximately  the  same,  i.e.,  .007  rads/ 
sec  for  both  situations.   However,  as  might  be  expected  from 

a  human,  the  turn  rate  exhibits  a  much  more  sinusoidal  be- 
havior as  compared  to  the  optimal  solution  which  immediately 
sought  and  settled  down  to  a  steady  state  value. 

3.  The  state  variables  under  real  time  simulation  showed 
a  more  pronounced  fluctuation. 

4.  An  advantage  of  the  real  time  simulation  over  the 
optimal  solution  was  that  the  human  operator  in  all  cases, 
i.e.  both  thrust  failure  and  thrust  reversal,  was  able  to 
maneuver  the  ship  back  to  the  original  track,  whereas  the 
optimal  solution  maneuvered  the  ship  to  a  parallel  track  and 
never  returned  to  the  original.   This,  of  course,  stems  from 
the  method  used  to  find  the  "steady  state"  values  around 
which  the  equations  of  the  system  were  linearized. 

Nevertheless,  the  difference  of  18  feet  between  the  original 
track  and  the  one  that  is  reached  at  steady  state  by  the  opti- 
mal algorithm  is  not  significant  compared  either  to  the  ship's 
width  (100  feet)  or  to  the  tactical  considerations  of  the  ship. 
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D.   NONLINEAR  VERSUS  LINEAR  SYSTEM  RESULTS 

In  order  to  determine  how  "optimal"  the  response  of  the 
nonlinear  system  was  with  the  use  of  the  state  feedback 
values  obtained  from  the  linear  regulator  design,  a  com- 
parison of  values  of  the  states  and  the  performance  indices 
were  made  during  a  failure. 

The  results  are  shown  in  Table  III.   As  can  be  easily 
verified,  the  differences  in  steady  state  values  obtained 
from  the  linearized  system  and  the  nonlinear  system  are 
quite  small,  ranging  from  .77%  to  1.77%,  with  the  exception 
of  the  steady  state  thruster  angle  (15%  difference) . 

The  above  results  point  out  that  the  linear  assumptions 
and  approximations  used  in  this  study  hold  quite  well  over 
the  operating  range. 

TABLE  IV 
Comparison  of  Linear  and  Nonlinear  Models 

LINEAR   NONLINEAR 
Surge  velocity  (ft/s)  75.98    74.63 

Sway  velocity  (ft/s)  h. 88     1.87 

-5  -5 

Angular  velocity  (rad/s)  .214x10     -.211x10 

Yaw  angle  (rad)  -.0248  -.025 

Displacement  in  x  direction  (ft)  7453  7342 

Displacement  in  y  direction  (ft)  10.762  10.88 

Steady  state  thruster  angle  (rad)  ,15256  .1804 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

Comparison  of  the  results  of  all  three  methods  examined 
leads  to  the  conclusion  that  the  optimal  control  solution  is 
by  far  the  best  one. 

Control  effort,  be  it  thruster  angle  or  thruster  magni- 
tude variation,  is  less  with  the  optimal  control  solution. 
This  is  of  great  consideration  not  only  from  the  aspect  of 
wear  and  tear  of  the  mechanical  parts  but  also  from  the 
point  of  ship  performance  as  a  naval  unit,  especially  since 
its  speed  is  so  far  greater  than  the  conventional  ships . 

Augmentation  of  human  operator  control  inputs  with  the 
optimal  solution  during  failure  scenarios  presents  the  best 
solution  in  order  to  maintain  positive  control  of  the  ship 
under  any  operating  conditions.   This  approach  would  combine 
the  best  attributes  of  the  human  operator  and  the  optimal 
control  schemes. 

B.  RECOMMENDATIONS 

Based  on  the  previous  analysis,  the  following  could  be 
investigated  in  the  future. 

1.   A  sufficiently  complex  constrained  parameter  optimi- 
zation routine  (COPES-CONMIN)  can  be  used  to  arrive  logically 
at  a  set  of  relative  system  response  with  respect  to  physical 
constraints.   This  optimization  will  still  be  dependent  on  the 
performance  index  selected  by  the  design  engineer  and  his 
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his  ability  to  calculate  explicitly  its  value  for  a  given 
dynamic  response  [Ref.  11]. 

2.  The  right-hand  side  of  the  equations  of  motion  could 
be  expanded  to  include  more  terms  to  provide  an  increase  in 
quantitative  accuracy  of  the  drag  terms  such  as  bubble  drag 
and  include  also  a  constant  power  propulsion  plant  versus  a 
constant  thrust  one.   This  could  be  done  easily  since  the 
terms  are  additive. 

3 .  Development  of  a  computer  program  to  augment  the 
control  efforts  commanded  by  a  human  operator  should  be 
investigated,  tested  and  implemented  on  the  RTS5D  of  Ref.  4. 

4- .   It  might  prove  worthwhile  to  investigate  more  closely 
the  trial  and  error  design  for  improved  response  character- 
istics of  this  method  of  control. 
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APPENDIX  A 
LINEAR  REGULATOR  DESIGN  THEORY 

Optimal  control  theory  provides  a  means  for  utilization 
of  state  variable  feedback  in  order  to  control  a  system  by 
minimizing  a  performance  index  or  cost  function  [Ref s .  7  and 
8]. 

Having  a  linear  system  defined  as  x  =  Ax  +  By   where  y 
is  the  control  input  vector,  the  linear  op  Imal  control  law 

is  defined  as  u   .(t)  =  -Gx(t) . 

~opt       ~~ 

The  assumptions  under  which  this  solution  is  valid  are: 

1.  The  control  vector  is  unconstrained. 

2.  The  performance  index  is  defined  as 

J  =  /tf   {xTQx  +  yTRy>  dt 
o 

where  Q  is  a  positive  semidefinite  or  positive  definite  real 
symmetric  matrix  and  R  is  a  positive  definite  real  symmetric 
matrix.   The  matrices  Q  and  R  determine  the  relative  impor- 
tance of  the  square  of  the  error  and  energy  expenditure 
during  the  transient  period. 

3.  The  system  is  controllable  and  state  observable. 
With  these  assumptions,  substitution  for  y  leads  to 

J  =  /  f   (xTQx  +  G(t)TxT  R  G(t)x)  dt 


J  =  /tf{xT  (Q  +  G(t)T  R  G(t))x}  dt 


84 


The  matrix  G(t)  is  found  to  be 

G(t)  =  R"1BTP(t) 
Here  P(t)  satisfies  the  matrix  Riccati  equation 

P(t)  +  Q  -  P(t)  BR"1BTP(t)  +  P(t)A  +  ATP(t)  =  0 


subject  to  the  boundary  condition  P(tf)  =  0. 

The  elements  of  the  matrix  G(t)  are  referred  to  as  the 
feedback  coefficients,  since  the  optimal  control  consists  of 
a  time-weighted  linear  combination  of  the  state  variables. 

This  result  is  important  since  it  indicates  that  in  an 
optimal  system  all  the  state  variables  are  to  be  fed  back, 
not  only  the  output,  as  it  is  customary  in  classical  control 
theory.   A  graphical  representation  of  this  optimal  system 
is  presented  in  Fig.  42. 

The  Riccati  equation  is  integrated  backwards  in  time  from 
the  known  terminal  condition  P(  °°  )  =  0  until  the  desired 
degree  of  accuracy  is  reached.   When  the  final  time  approaches 
infinity,  the  P  matrix  becomes  constant,  thereby  simplifying 
its  determination. 

As  can  be  seen  for  the  previous  presentation,  a  different 
linear  optimal  control  is  generated  for  every  relative  set 

of  weighting  matrices,  Q  and  R,  provided.   Additionally,  the 

pi      <> 

system  response  is  optimal  in  a  mathematical  sense  only,  since 
no  assurance  on  acceptable  system  overshoots  and  input  peak 
values  can  be  guaranteed. 

A  minimum  area  could  have  an  excessive  peak  which  would 
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provide  a  physically  unrealizable  control  scheme.   The 
design  synthesis  thus  requires  an  initial  guess  at  the  rela- 
tive importance  multipliers,  followed  by  subsequent  adjust- 
ments based  on  time  simulation  of  the  controlled  system 
dynamics  . 
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